Submitted to the Astrophysical Journal. 



A Viewing Angle - Kinetic Luminosity Unification 
Scheme For BL Lacertae Objects. 

Markos Georganopoulos and Alan P. Marscher 
Department of Astronomy , Boston University, 725 Commonwealth Avenue, Boston, MA 

02215 

ABSTRACT 

We propose a unified classification for BL Lac objects (BLs), focusing on 
the synchrotron peak frequency u s of the spectral energy distribution. The 
unification scheme is based on the angle G that describes the orientation of the 
relativistic jet and on the electron kinetic luminosity Kkin of the jet. We assume 
that Akin scales with the size of the jet r in a self-similar fashion (Ak in oc r 2 ), 
as supported by observational data. The jets are self-similar in geometry and 
have the same pressure and median magnetic field at the inlet, independent 
of size. The self-similarity is broken for the highest energy electrons, which 
radiate mainly at high frequencies, since for large sources they suffer more severe 
radiative energy losses over a given fraction of the jet length. We calculate 
the optically thin synchrotron spectrum using an accelerating inner jet model 
based on simple relativistic gas dynamics and show that it can fit the observed 
infrared to X-ray spectrum of PKS 2155-304. We couple the accelerating jet 
model to the unification scheme and compare the results to complete samples of 
BLs. The negative apparent evolution of X-ray selected BLs is explained as a 
result of positive evolution of the jet electron kinetic luminosity Afc in . We review 
observational arguments in favor of the existence of scaled-down accretion disks 
and broad emission-line regions in BLs. The proposed unification scheme can 
explain the lack of observed broad emission lines in X-ray selected BLs, as well 
as the existence of those lines preferentially in luminous radio-selected BLs. 
Finally, we review observational arguments that suggest the extension of this 
unification scheme to all blazars. 



Subject headings: radiation mechanisms: non-thermal — BL Lacertae objects: 
general — galaxies: jets 
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1. INTRODUCTION 

The family of blazars includes those AGN (active galactic nuclei) that are characterized 
by compact radio morphology and variable, polarized nonthermal continuum. Apparent 
superluminal motion and 7-ray emission are also common properties among blazars (for 
a review of their properties see Urry & Padovani 1995). When the spectrum exhibits the 
usual quasar-like broad emission lines, the object is classified as a flat spectrum radio 
quasar (FSRQ). On the other hand, the term BL Lacertae object (BL) is reserved for the 
lineless or almost lineless blazars. The parent population of FSRQs ( Padovani 1992|) and 
BLs ( [Padovani fc Urry 1990| ) is thought to consist of Fanaroff-Riley type II (FR II) and 
Fanaroff- Riley type I (FR I) radio galaxies ( [Fanarotf fc Riley 1974 ), respectively. 



The standard interpretation of the nonthermal continuum radiation of blazars from 
radio to 7 rays is synchrotron and inverse Compton emission from a collimated relativistic 
plasma jet ( piandford fc Rees 1978| ). The currently popular paradigm involves the existence 
and dynamical consequences of a magnetic field that threads an accretion disk around a 
massive black hole ( |Blandford fc Payne 1982j ). The collimation and acceleration of the jet 
in such a scenario is currently a field of active research; both analytical (e.g. [Appl 1996]) 



and numerical (e.g. Puyed, Pudritz, fc Stone 1997|) work shows promising results thus far. 



The observational taxonomy of BLs consists, according to the method of discovery, 
of radio selected (RBLs) and X-ray selected (XBLs) sources. According to Padovani & 
Giommi (1995), this corresponds to a more physical dichotomy, based on the frequency at 
which the peak of the synchrotron spectral energy distribution (SED, quantified by vL u ) 
occurs: low-frequency peaked BLs (LBLs), which are mostly RBLs, and high-frequency 
peaked BL Lac objects (HBLs), which are mostly XBLs. The XBLs are less variable and 
polarized than RBLs flJannuzi, Smith, fc Elston 1994] ) and have flatter optical spectra 



( [Ledden fc 1 Dell 1985| ). The RBLs are more core-dominated flPerlman fc Stocke 1993| ) 



and for a given X-ray luminosity are more luminous at radio and optical frequencies than 
XBLs ( pVIaraschi et al. 1986|) . RBLs are characterized by a weak positive evolution (i.e., 
dependence of number counts on redshift; |Stickel et al. 199 TP , while XBLs exhibit a negative 



evolution (|Perlman et al. 1996| ). The recent discovery of a population of intermediate BLs 



( Laurent-Muehleisen 1997| |Perlman et al. 1998|) suggests a continuous family of objects 



rather than two separate classes. 

Maraschi et al. (1986) proposed that the X-ray emission comes from a compact, less 
beamed region, while the radio emission comes from an extended, more beamed region. 
According to this "orientation hypothesis", the RBLs form a small angle (O <^ 10°) between 
the line of sight and the direction of motion of the emitting plasma, while the XBLs form a 
larger angle (10° <^ Q <^ 30°). A more recent interpretation of the differences between XBLs 
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and RBLs ( Padovani fc Giommi 1995|) is based on the cutoff frequency of the synchrotron 



SED. According to this "SED-cutoff hypothesis," most of the BLs are characterized by a 
cutoff in the IR/optical band, and these are the radio selected objects. The small fraction 
of BLs that have a cutoff at UV/X-ray energies are the X-ray selected BL Lac objects. This 
model does not offer a physical explanation for the different cutoff frequencies of the SED. 

There is an ongoing debate in the literature about the actual reason for the differences 
between XBLs and RBLs. Padovani Sz Giommi (1995) argue that the observed range of the 
core dominance parameter R (the ratio of core to extended radio flux) of the XBLs does 
not agree with the orientation hypothesis. On the other hand, Kollgaard et al. (1996a) 
conclude that the fact that the mean core radio power of RBLs exceeds that of XBLs by 
more than an order of magnitude for sources selected according to extended radio power, 
is incompatible with the SED-cutoff hypothesis. Comparison of VLBI images of RBLs 
and XBLs QKollgaard, Gabuzda. & Fcigclson 1996|) suggests that the jets in XBLs fade 



faster than in RBLs, which supports the idea that factors other than orientation are also 
important in explaining the differences between XBLs and RBLs. 

Urry & Padovani (1995) note that the SED-cutoff hypothesis cannot explain the 
different polarization properties of XBLs and RBLs, particularly the stability of the 
polarization angle in the XBLs. Lamer, Brunner, & Staubert (1996) analyze a sample of BL 
Lac objects observed with the PSPC detector on board the ROSAT satellite and note that 
the SED-cutoff hypothesis can explain the range of the crossover frequency between the 
soft and the hard components in the X-ray regime better than the orientation hypothesis. 
Sambruna, Maraschi, & Urry (1996), using complete samples of BLs, reach the conclusion 
that the range of the peak frequency of the synchrotron energy distribution in BLs cannot 
be reproduced under the orientation hypothesis. An alternative scenario, which connects 
the synchrotron peak luminosity L s to the synchrotron peak frequency u s of the SED of 
the BLs through an empirical relation, has been proposed recently by Fossati et al. (1997). 
One feature that this scenario reproduces more successfully than the previous two is the 
redshift distribution of XBLs and RBLs. 

In this work we propose that two parameters determine the observed characteristics 
of a BL. The first which must play a role, given the relativistic nature of the jet flow, is 
the angle O formed between the line of sight and the jet axis. The second is the electron 
kinetic luminosity A kin of the jet. We start by developing and testing a numerical code for 
the accelerating inner jet model, which we use as our working hypothesis for the physical 
description of the jet. Invoking recent observational studies, we propose a new classification 
scheme, whose main observational parameter is the synchrotron peak frequency u s of the 
SED of BLs. We then introduce the theoretical O-A unification for BLs, based on two 
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parameters: the angle G between the jet axis and the line of sight, and the electron kinetic 
luminosity A kin of the jet. These are coupled to a simple self-similar jet description that 
relates A kin to the size of the jet (A kin oc r 2 ). We use the accelerating jet model to show- 
how this unification scheme can be used to explain several characteristics of BLs, including 
the negative apparent evolution of the XBLs. We also show, using the self-similarity 
scenario we propose, that a picture in which BLs have a scaled-down accretion disk and 
broad emission-line region (BELR) environment similar to those observed in FSRQs and 
high polarization quasars (HPQs), is in agreement with observations. (We will use the 
terms FSRQ and HPQ interchangeably to signify all the non-BL blazars). Finally, we 
review observational evidence that supports a unified picture for all blazars under the above 
self-similar scheme. 



2. THE JET MODEL 

A first approach in calculating the synchrotron emission produced by an accelerating 
and collimating plasma flow was presented by Marscher (1980): In the accelerating 
inner jet model (based on the work of [Blandford fc Rees 1974| and [Reynolds 1982] ), 
continuously generated ultrarelativistic plasma is confined by pressure (hydrostatic and/or 
magnetohydrodynamic) that decreases along the jet axis. The internal energy of the plasma 
is converted to bulk kinetic energy and the jet is accelerated and focused. The electrons 
interact with a predominately random magnetic field and cool through synchrotron radiation 
and adiabatic expansion. At the same time, inverse Compton losses become important 
when the synchrotron photon energy density becomes comparable to the magnetic field 
energy density. Higher frequency radiation is emitted close to the base of the jet, since only 
there are the electron energies high enough to do so. Lower frequency photons are emitted 
there as well as farther downstream. The velocity of the jet increases with distance; this 
implies larger Doppler boosting for greater distances down the jet and therefore at lower 
frequencies, out to the point where the Lorentz factor T <^ B _1 . The relevance of this 
approach to the BL Lac phenomenon was strengthened through studies that suggested that 
the X-ray emitting region is less relativistically boosted than that of the radio emission 
( [Maraschi et al. 1986] ; |Urry, Padovani, fc Stickel 1991Q . 



We use the following phenomenological description for the accelerating inner jet: 
ultrarelativistic plasma (mean ratio of total to rest mass energy per particle 7 3> 1 in the 
proper frame of the fluid) is continuously injected at the base of the jet. The bulk flow 
of the plasma at the injection point is parameterized by the bulk Lorentz factor T*. The 
axially symmetric pressure profile, with the pressure on the symmetry axis being less than 
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the equatorial pressure, defines a preferred direction for the expansion of the plasma. The 
plasma flow is progressively accelerated and collimated, converting its internal energy to 
bulk kinetic energy. The adopted phenomenological description of ultrarelativistic particle 
injection at the base of the jet can be linked, for example, to particle acceleration in a 
standing shock front in the plasma flow, formed due to a drop in the confining pressure 
( Pomez et al. 1997|) . Since the emission at higher frequencies is confined to a region closer 
to the base of the jet than at lower frequencies, we expect ( [Marscher 1980| ; |Ghiselliiil 
fc Maraschi 19891) shorter variability timescales at higher frequencies, given that the 
acceleration of the flow is mild. In addition, the Doppler effect is comparatively more 
important for the lower than for the higher frequencies. Therefore the jet orientation affects 
the lower frequencies more than the higher ones. 



2.1. Flow description 

We adopt a flow description (see appendix |A|) similar to that of Blandford & Rees 
(1974) and Marscher (1980), the latter of which considered as the base of the jet the sonic 
point, where the bulk Lorentz factor has the value = \JZJ2. If we identify the base of 
the jet instead with a standing shock, we can relax the assumption that the flow is sonic 
at the base of the jet, and allow for higher values of the bulk Lorentz factor I\ at that 
point. The minimum distance of the base of the jet from the accretion disk is obtained 
by requiring that the electron Thomson losses due to the accretion disk photons are less 
significant than the synchrotron losses (see Appendix |B|). The quantities that describe the 
flow are the total (electron + proton) kinetic luminosity At of the jet, the radius r* of the 
base of the jet, and the exponent e and the length scale z± that describe how fast the jet 
opens and accelerates. Assuming energy equipartition between electrons and protons, the 
electron kinetic luminosity Akin is 

A km = A t /2. (1) 

Thus half of the injected energy is in protons that do not radiate, and therefore we can 
reasonably approximate the flow as adiabatic. 

The bulk Lorentz factor as a function of distance along the jet axis is given by equation 

r W = r.(£)\ (2 ) 

while the radius of the jet is given by a modified form of equation (|A9|) : 

i 

4 

(3) 
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For 3-D steady flow, this equation applies to each streamline, but with r(z) replaced by r str . 
The flow velocity is tangent to the streamlines and the loci of constant scalar quantities such 
as density are surfaces that are perpendicular to the streamlines. Note that for I\ = \Js/2 
equation (|3|) is reduced to equation 



2.2. Magnetic field 



For the purpose of this work, we will consider the magnetic field to be mainly random 
in direction, with a small ordered component parallel to the jet axis. Under the assumption 
of an isotropic electron distribution (see §§ |2.4| ), the choice of magnetic field geometry is not 
critical, since we are interested here in the total and not the polarized flux. Using simple 
flux conservation arguments, we assume that the magnetic field in the fluid frame decays as 
B oc r~ m , 1 < m < 2, or in terms of the jet axial coordinate z: 



Biz) 



11 ( - 

.2* 



(r2 - 1; 



2c 



(4) 



2.3. Electron energy losses 

We assume that the electrons are injected at the base of the jet and propagate 
downstream while losing energy due to adiabatic and synchrotron energy losses, without 
any reacceleration (except at downstream shocks fronts, the radiation from which is not 
considered here) . The electrons will also suffer inverse Compton losses from the synchrotron 
photons they produce (synchrotron self Compton (SSC) model, e.g., [Marscher fc Bloom| 
1994 ) as well as from photons coming either from the broad-line clouds flSikora, Bcgclma'n] 



|fc Rees 1994| ) or from the putative accretion disk (Permer fc Schlickeiser 1993|) (external 
Compton (EC)). In our treatment we require that the base of the jet is sufficiently far from 
the central engine, so that the electron Thomson losses due to scattering of the accretion 
disk radiation are smaller than the synchrotron losses (see Appendix B). A self-consistent 
treatment of inverse Compton scattering is extremely complicated, since the local photon 
energy density results from photons that have been emitted in different regions of the jet 
at different retarded times. Observations using the Compton Gamma Ray Observatory 
(CGRO) ( |von Montigny et al. 1991% Pondi fc Ghisellini 1995| ; |Fossati et al. 199fl ) suggest 



that during flares the Compton luminosity is usually greater than the synchrotron in 
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FSRQs by rs 1-2 orders of magnitude, but is comparable to, or less than, the synchrotron 
luminosity in BLs, and during low states of 7-ray activity the synchrotron luminosity 
dominates the power output in BLs. In SSC models the beaming cone of the synchrotron 
and Compton radiation is the same, and the ratio of apparent luminosities of the two 
components equals the ratio of synchrotron to Compton electron losses. 

In EC models the beaming cone of the inverse Compton radiation is more narrow 
than the synchrotron cone ( |Dermer 1995| ). Statistically, in the case of EC models, even if 



Compton losses are equal to the synchrotron losses, in any BL sample selected at frequencies 
dominated by synchrotron emission, one would expect that the Compton luminosity would 
be significantly greater than the synchrotron one for a fraction of the sources equal to the 
ratio of the beaming cone solid angles. Since this is not observed, it seems plausible to 
assume that the Compton losses even in EC dominated BL models are not the main energy 
loss mechanism for BLs. 

We combine the adiabatic and synchrotron losses to obtain the electron energy-loss 
equation in the fluid frame: 



where 



dz 



'dnfjz)' 

dz 



-erf 



+ 



synch 



dz 



ad 



-3me 



2c 



n - -1 



X 



2 p 4 



771—1 

2 



£7 

z 



(5) 



(6) 



Since (dj / dz) sync h oc 7 2 while (d'y/dz) a d oc 7, synchrotron losses dominate for the 
higher energy electrons. The Lorentz factor 7 describes the internal kinetic energy of the 
electrons in the fluid proper frame. Solving for the electron energy, we obtain 



where the integral 



7(z) 
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z / z \ -(3m+l)e 
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must be solved numerically for m^l, 



2.4. Electron energy distribution 

The EED is, in principle, determined by the particle acceleration and the dominant 
energy loss mechanisms acting simultaneously on the electrons ( |Kirk 1997]) . These two 
factors should determine the shape and the cutoffs of the distribution. In our model we 
assume that the initially injected EED has a power law form: 

iV^,^) = N^~ s , 7* G [7min(^f),7mox(«*)]- (9) 

We further assume that the initial electron energy distribution (EED) remains isotropic 
in the comoving frame, through redistribution of the electron pitch angles due to their 
scattering on plasma waves. The EED distribution is therefore a function of energy and 
position along the jet axis only. The normalization factor is related to the comoving 
energy density e e /* at the base of the jet 

eeU = / 7^7 "dl = TT—iTmal ~ Train)- (10) 



Using equation ( |A2| ) for the base of the jet (A = r 2 .) and the relation e* = (Z/4)w± (see 



Appendix |A|), together with the assumed energy equipartition between electrons and 
protons (e e z* = e*/2), we obtain 

N = 3 ( 2 ~ g ) Afc - (in 

4T 2 8+ch 2 - s - 7^ SN ) 1 J 

Imax IrninJ' 

In our treatment it is implicitly assumed that the acceleration time is much shorter 
than the radiative loss time for the considered frequencies, and we can therefore separate 
the region where acceleration takes place from the radiative zone. This assumption may 
break down for the highest synchrotron frequencies, although even in extreme objects like 
Mkn 421 the 1-10 keV variability behavior indicates that this assumption is valid QKirkj 



1997|) . We therefore consider this approach adequate for synchrotron radiation with photon 
energy up to ~ 10 keV. Conservation of particles and eq. (0) give the evolved energy 
distribution: 



N( 7 ,z) = N^- s (±) 



z \ -e(s+2) 

X 



z x ■ X 8-2 



1-7(^)^(^-1)-^)) . (12) 
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This profile is no longer a power law, and on a log-log plot it gradually steepens toward 
higher frequencies. Note also that the EED evolves as the plasma moves downstream. This 
gives rise to a non-power law locally emitted synchrotron spectrum in the local flow frame, 
which changes with z in frequency cutoff and shape. 



2.5. Relativistic transformations 

Given the magnetic field and the EED in the local fluid proper frame we can calculate 
the synchrotron emission coefficient in the fluid frame using standard synchrotron formulae 
flPacholczyk W?§ : 

x / N( 1 ,z)F(x)d 7 , (13) 

where 

poo 

F(x) = x / K 5/3 (z)dz. (14) 

J X 

K 5 / 3 is the modified Bessel function of order 5/3 and x a dimensionless frequency: 

v 

x 



where 

3e 

Vcrit(z) = B(z) smcj)'^(z) 2 . 

In eq. ([H|) is the angle between the magnetic field and the line of sight in the fluid frame. 
In order to calculate this angle we first need to transform the angle 9 between the line of 
sight and the local flow velocity from the observer's frame to the fluid frame using standard 
relativistic aberration formulae ([Rybicki fc Lightman 1979]) : 



cos#-/3 
1 — p cos 

Having now the angle 9' between the line of sight and the local flow velocity in the 
local fluid frame, we can obtain 0' from the angle between the magnetic field and the local 
flow velocity in the fluid proper frame. 
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To transform the emission coefficient to the observer's (unprimed) frame we need to 
take into account the transformations between this frame and the fluid (primed) frame. If 
z is the redshift and 5 = 1/(7(1 — (3 cos 9)) is the Doppler boosting factor, where 9 is the 
angle between the sight and the fluid velocity, then the frequency in the observer's frame is 
v = u'5/(l + z) and the emission coefficient is ([Rybicki fc Lightman 19791) : 

^(iTl)^" (16) 

Having now the synchrotron emission coefficient in the observer's frame, we can perform 
the radiative transfer and calculate the total synchrotron flux for optically thin frequencies. 



3. MODEL SPECTRAL ENERGY DISTRIBUTIONS 

3.1. The parameter space 

We have performed an extensive study of the model parameter space in order to 
examine the response of the resultant SED to a variation of each physical parameter about 
the fiducial values given in Table [l]. The meaning of the symbols is the same as in [Q. In 
Figure [1] we allow one of the model parameters to vary in each panel and plot the resulting 
SED. 

As e increases (corresponding to a more poorly collimated jet) the peak apparent 
synchrotron luminosity L s decreases. This is due to an increase of the adiabatic losses and 
a decrease of the synchrotron power (cf. eq. (5|), regardless of the value of m. At the same 
time the peak frequency v s of the SED increases, since the synchrotron losses decrease, and 
the effective break of the EED shifts to higher energies. An increase in the magnetic field 
amounts to an increase in L s and a decrease in v s . This is expected since a higher magnetic 
field results in higher synchrotron power and hence losses, primarily for the higher energy 
electrons. 

The behavior of the SED as we increase is less straightforward. As I\ increases, 
v a increases, the IR luminosity decreases, the X-ray luminosity increases, and L s increases 
until I\ exceeds a critical value, beyond which it decreases. The key in understanding this 
behavior is that the jet emission is stratified, with the higher frequencies, which emerge 
closer to the base of the jet, being characterized by smaller Lorentz factors. As we increase 
I\, the effective Lorentz factor T u for each frequency increases, and the angle (~ 1/T„) 
into which the radiation is beamed decreases. In the IR regime (for the parameters chosen 
in Table |TJ) an increase in I\ beams most of the radiation into a cone with an opening 
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angle smaller than O; this causes the Doppler de- amplification we observe. In X-rays the 
increase in I\ beams the radiation into a narrower cone that still includes the line of sight 
and therefore the observed luminosity is more Doppler boosted. This frequency-dependent 
response causes u s to shift toward higher frequencies. L s increases as long as the Doppler 
boosting of the higher frequencies can compensate for the Doppler de-amplification at the 
lower frequencies. After a certain point (for the parameter values we use, this occurs at 
I\ « 3), the observed peak luminosity declines, since most of the emitted power is beamed 
into a cone that does not include our line of sight. 

An increase in z* produces a behavior similar to that exhibited by the increase in B*. 
This is expected, since an increase in z* causes the magnetic field to decline more slowly 
with time in the electron's frame. That is, the electron spends more time in the roughly 
uniform magnetic field region as the size of that region increases. 

An increase in s (a steeper EED) produces a steeper spectrum, as expected. As the 
EED becomes steeper, the energy content of the high energy tail of the EED decreases. 
Since these are the electrons that lose energy faster and emit at higher frequencies, both L s 
and u s decrease as the EED becomes steeper. An increase in the upper cutoff ^ ma x of the 
EED affects dramatically the high frequency tail of the SED, making the X-ray spectrum 
flatter and increasing both L s and u s . The effect on the low frequencies is minimal, as 
expected, since the high energy electrons radiate a small fraction of their power there. The 
behavior of the SED resembles that described by Padovani & Giommi (1995), who attribute 
the cause of the differences between XBLs and RBLs to different values of 7 ma x- 

The exponent m, which describes the evolution of the comoving magnetic field along 
the jet axis, and the radius of the jet (not shown in Figure only weakly affect the 
SED. Increasing m from m=l to m=2 decreases the IR luminosity by a factor of ~ 2, 
while the X-ray luminosity remains essentially the same. Increasing r* by a factor of four 
reduces the X-ray luminosity by a factor of ~ 2, while the IR luminosity remains constant. 
Note, however, that in a time dependent simulation the minimum variability timescale 
[t V ar > r*sinO/c] provides an upper limit to r*. The value of the lower cutoff 7 min of 
the EED should in principle be derived from the particle acceleration theory. For a given 
jet electron kinetic luminosity, as 7 m i„ increases, the particle density and the observed 
synchrotron luminosity increase, since the average energy per particle increases. Hence, Akin 
and 'jmin are closely connected and only one (A kin ) needs to be treated as a free parameter 
unless the value of the density is needed. For the purpose of this work we will assign the 
values m—l and 7™ n = 10 2 . 

For a given physical description of the jet, the angle between the jet axis and the 
line of sight is still a free parameter, in the sense that the jet can be oriented at any 
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angle relative to the line of sight. As mentioned in the introduction, one of the leading 
scenarios for the interpretation of the differences between RBLs and XBLs is based on the 
orientation of the jet relative to the line of sight. We have run our simulation using the 
parameters of Table |I], and have allowed the angle to vary, as presented in Figure ^. The 
data points refer to the 1 Jy sample of RBLs (triangles) and the Einstein Extended Medium 
Sensitivity Survey sample of XBLs used by Sambruna et al. (1996). As can be seen in the 
bottom panel, as the angle increases from = 0° to = 20° the synchrotron luminosity 
L s decreases by a factor of ~ 20 and the synchrotron peak frequency v s increases by ~ 2 
orders of magnitude. At the same time the spectrum becomes flatter at all frequencies. 
The synchrotron peak luminosity L s (middle panel) and the optical to X-ray spectral index 
a ox (top panel), calculated taking into account the luminosities at 5500 A and at 1 keV, 
are plotted versus the peak frequency u s of the SED. Note that the model point shifts 
from a predominately RBL region to a predominately XBL region in both the luminosity 
and spectral index plots, as the angle increases. This behavior is not specific to the set of 
model parameters used here, but appears in a qualitatively similar manner for a wide range 
of initial conditions: As the angle increases, L s decreases, v s increases, and the spectrum 
becomes flatter. Although these differences are similar to the differences found between 
XBLs and RBLs, previous work (Sambruna et al. 1996; Georganopoulos fc Marscher 1996[ ) 
has shown that a change only in the jet orientation is not enough to explain the range of 
observed properties of BLs, hence additional physical differences must be invoked. 



3.2. Simulation of the SED of PKS 2155-304 

As an application of the model we simulate the SED of PKS 2155-304, a relatively 
nearby ( redshift .2=0.116) bright XBL. Its 7-ray luminosity is less than its synchrotron 
( |Vestrand, Stacy, fc Sreekumar 1995| ), which implies (see also section |2.3| ) that synchrotron 



losses are more important than inverse Compton losses. The object has been the target of 
two intensive multiwavelength campaigns, one during November 1991 ( Edelson et al. 1995| ). 



and one during May 1994 (|Urry et al. 1997] ). 



Since PKS 2155-304 is highly variable and a strictly simultaneous SED is practically 
unattainable, we must approximate a simultaneous SED with one constructed from fluxes 
measured during a time interval smaller or at least comparable to the timescale of the 
fastest significant variations. The shortest variability timescale for the November 1991 
dataset was ~ 0.7 days. We construct the SED plotted in Figure [3| using observations in the 
time interval between November 14.0 and 14.76. This time interval is approximately equal 
to the fastest variability timescale observed, and it corresponds to a local flux minimum for 



13 



the optical to X-ray frequency range. The solid curve in Figure |3| represents the best-fit 
model. The values of the parameters used are given in Table |2|. The \ 2 f° r this fit is 49.65 
for nine degrees of freedom. This high value can be partly attributed to the arbitrariness 
of the selected observed fluxes used to construct the "base" observed SED, but in any case 
is typical of fits to multiwavelength spectra with non-complex models. Since an exhaustive 
exploration of the parameter space is practically impossible, we cannot assert that this is 
the best possible fit. To derive a more constrained description of the source, we would need 
to model successfully the observed variability behavior. 



4. A NEW SED-BASED CLASSIFICATION SCHEME 

BL Lac objects are discovered mostly through X-ray or radio surveys. The existing 
complete samples therefore are either X-ray ( [Morris et al. 199T| |Perlman et al. 1996|) or 



radio flux-limited samples ( (Stickel et al. l99~T|) and BLs are named XBL (X-ray selected) 
or RBL (radio selected), depending on the method of discovery. In general RBLs have 
steeper spectra than do XBLs. The SEDs of RBLs peak somewhere in the infrared and 
have a higher synchrotron peak luminosity L s than do the XBLs, which usually peak in the 
UV-soft X-ray range (Sambruna et al. 1996). 

A recently proposed classification scheme for BLs (Padovani, Giommi, fc Fiore 1996] ) 
corresponds more closely to the observed characteristics of the objects than the method 
of discovery. According to this scheme BLs are classified as high frequency peak BLs 
(HBL) or low frequency peak BLs (LBL), depending on the value of the X-ray to radio flux 
ratio. Objects with f x /f r < 10~ 1L5 are classified as LBLs and objects with f x /f r > 10 -1L5 
as HBLs. Most of the HBLs are XBLs and most of the LBLs are RBLs. Although this 
classification scheme is based on observed properties of BLs, it inherits the dichotomization 
of BL Lac objects that was imposed by the two different discovery methods. Such a 
dichotomization can hinder our efforts to understand the physics of BLs, particularly if the 
two methods of discovery force us to observe the two extremes of a continuous distribution, 
thereby introducing a selection-induced bimodality. 

Strong indications that a bimodality is introduced by the observational techniques 
come from a recent study by Laurent- Muehleisen (1997). The cross-correlation of the Rosat 
All Sky Survey (RASS) and the Green Bank 5 GHz radio survey revealed the existence 
of a previously highly undersampled population of BLs (120 objects) with properties 
intermediate between those of XBLs and RBLs. Specifically, the X-ray to radio flux ratio 
of these sources and their location on a diagram of radio to optical (a ro ) versus optical 
to X-ray (a ox ) spectral index bridge the gap between XBLs and RBLs. Similar results 
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( [Pcrlman et al. 1998 ) are emerging from the deep X-ray radio blazar survey (DXRBS). 



The question that naturally arises is whether there is an observed quantity that varies in a 
continuous fashion, has a well defined physical significance, and can be used as a classifying 
parameter. We argue that such a quantity is the synchrotron peak frequency v s . 

Recently, Sambruna et al. (1996) examined the multifrequency spectral properties 
of three complete samples of blazars: the Einstein Extended Medium Sensitivity Survey 
sample (EMSS) of XBLs (23 sources; |Morris et al. 1991| ), the 1 Jy sample of RBLs (29 



sources; |Stickel et al. 1991| ), and a small complete sample of FSRQ from the S5 survey 
(8 sources; [Brunner et al. 1994| ). After they calculated the rest-frame fluxes and spectral 



indices for each object, they fit a parabola to the synchrotron power {vL u ) spectrum. X-ray 
fluxes were excluded from the fitting in the cases in which the X-ray spectral index a x 
was flatter than the broadband optical to X-ray spectral index a ox , a condition that is 
thought to indicate that inverse Compton radiation dominates over the steep high-energy 
synchrotron tail of the spectrum. In Figure |]a-f we plot the spectral indices a ox , a ro , 
a rx , the X-ray concavity index a ox — a x , and the synchrotron bolometric Lb and peak L s 
luminosities as a function of the synchrotron peak frequency v s for the two BL samples of 
Sambruna et al. (1996). If we consider the two samples collectively, we find that in all six 
diagrams there is an upper envelope separating the populated area from a well-defined zone 
of avoidance. The upper envelope is such that for a given peak frequency v s there is a range 
of "avoided" luminosities and spectral indices. As v s increases, the maximum observed 
luminosity decreases, the steepest observed spectral index flattens, and the maximum 
observed concavity index a ox — a x decreases. On the other hand, the flattest spectral 
indices observed and the most negative concavity indices do not seem to be very sensitive 
to the value of v s . 

We propose a continuous classification system for BLs, based on the synchrotron 
peak frequency v s of the synchrotron SED. For example, a BL peaking at log(^ s ) ~ 14 
will be classified as a BL14. The class of a BL Lac object then defines the maximum 
luminosity of the object, the maximum steepness of its synchrotron spectrum over a range 
of frequencies, and the maximum concavity index. We note here that (unlike the transient 
variability mentioned before) long term, persistent changes in u s , should be reflected in the 
classification of the object. 

The synchrotron peak frequency v s of a BL is a very important parameter from a 
physical point of view as well. It is the frequency in which we receive the highest synchrotron 
power per logarithmic frequency interval. It is also closely related to the peak of the inverse 
Compton emission in both the SSC (Bloom fc Marscher 1996]) and the EC models ( (Bikora 



Begelman, fc Rees 1994j ; Permer fc Schlickeiser 19"93|) . This frequency is closely linked to 
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the the energetics and geometry of the synchrotron source and to the angle formed between 
the observer and the plasma velocity, if the source is moving relativistically (see §[|). 

As Sambruna et al. (1996) point out, the determination of v s requires simultaneous 
multiwavelength flux measurements and can be particularly sensitive to a lack of flux 
measurements close to the peak of the distribution. The problem of multiwavelength 
coverage becomes more difficult when v s happens to be either in the observationally 
inaccessible frequency range between sub-mm and infrared frequencies v w 10 12 — 10 14 Hz, 
or in the UV due to galactic absorption. If v s is in the X-ray energy range, and there are 
no flux measurements to constrain the parabola at frequencies higher than i/ s , the error in 
v s can be very large. This is the case for the two sources that have \o<gv s > 18, the XBL 
1433.5+6369, and the RBL 2005-489 (see Fig. 4 of Sambruna et al. (1996)), which we 
decided not to consider further in our study. In general, the error in the determination of 
v s should be smaller for objects with SEDs peaking at IR/optical frequencies, where v s is 
more easily constrained. 

We note, however, that even though v s is difficult to measure exactly, an approximate 
estimate within a decade of frequency is both feasible and adequate for the purpose 
of classification. Furthermore, the accuracy can be improved if one decides to perform 
multiwavelength observations with this specific goal in mind, instead of using the available 
data from the literature. 

Variability can also affect v s . In the case of individual flares, initially the spectrum 
hardens ( [Ulrich, Maraschi, fc Urry 1997] ) and the peak synchrotron frequency increases. 



This change is temporary, and the source returns to its pre-flare characteristics after the 
end of the flare. The most extreme recorded case of such a shift occurred in Mkn 501, that 
increased its synchrotron peak frequency by more than a factor of 100, changing form a 
BL17 to a BL19 flPian et al. 1998] ). Usually, though, the shift in v Sl something neither 
easily nor customarily calculated in multiwavelength studies, is approximately a factor of 
10 in the case of the few well documented, large amplitude flares. In the case of Mkn 421 
during the May 1994 (Macomb et al. 1995|) and the April-May 1995 QBuckley et al. 1996|) 



campaigns, v s increased by a factor « 10-100, while during the PKS 2155-304 May 1994 
campaign, v s increased by less than a factor of 10 ( |Urry et al. 1997|) . Since these excursions 



of v s are usually of the same magnitude as the uncertainty in measuring u s , and also of a 
transient character, we do not consider them a serious problem in assigning a specific v s to 
a BL, regardless of variability. 

Under our proposed unified classification based on v s , the questions that can direct us 
toward a physical unification are: (1) which are the physical and/or geometrical parameters 
that, given their range of values, create the range of the observed phenomenology? (2) Can 
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a model incorporate these parameters and derive the properties of complete samples of BL 
Lac objects? In the following section we propose a simple unification scheme based on a 
self-similar jet description. 



5. A 0-A UNIFICATION SCHEME 

5.1. A self— similarity hypothesis 

In order to the extended region that BLs occupy in the diagrams of Figure ^, we 
need to employ a model with at least two free parameters. Any model with only one free 
parameter can at best trace a curved trajectory in the observational graphs when this free 
parameter is allowed to vary. Under the assumption that the jets of BLs are relativistic, a 
change in the angle between the line of sight and the jet axis will induce a model-dependent 
change in the observed properties of the source. It therefore seems necessary to postulate 
minimally the existence of a physical free parameter that, together with the angle, is 
capable of reproducing the observed range of properties of BLs. 

The concept that relativistic jets can be described in a scale-invariant way using 
a self-similar description is very appealing, since these physical configurations span a 
luminosity range greater than ten orders of magnitude, from the jets in stellar systems in 
our galaxy to the most extreme FSRQs . Maraschi & Rovetti (1994), after studying samples 
of BLs, FSRQ, and FRI and FRII radio galaxies, suggested that all these sources operate 
in a self similar way, with the jet luminosity and angle between the line of sight and the 
jet axis accounting for the different observed properties. Gear (1993) studied the infrared 
colors of two complete samples of RBLs and XBLs, and hinted toward a BL unification 
based on the luminosity and the angle between the line of sight and the jet axis. 

We propose that the second, physical unification parameter is the electron kinetic 
luminosity of the jet Akin- At first view, one would be tempted to assume that A^n oc r*, 
since the Eddington luminosity scales with the black hole mass, which is proportional to the 
gravitational radius of the black hole. However, tying the luminosity of the system to the 
black hole mass cannot address the evolutionary characteristics of BLs. The black hole mass 
slowly increases, due to accretion of fresh material, and this translates into a slow increase 
of the luminosity, which essentially translates to a weak negative cosmological evolution, 
i.e., each source becomes more powerful as cosmic time increases. We know that this is the 
case neither for RBLs (see §§ [5.3|) , nor for all the other radio-loud AGN families (|Urry, fc 



Padovani 1995| ), which we presume are described by the same paradigm — a massive black 



hole, an accretion disc, and a relativistic jet. 
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An alternative description has been proposed by Rees et al. (1982). Under this 
scheme, the electromagnetic power L EM extracted from of a rapidly rotating black hole is 
proportional to a 2 , where a < r g is a length that measures the specific angular momentum 
ac of the black hole, and r g is the gravitational radius of the black hole. Assuming now that 
Lem and a are related linearly to A^n and r* correspondingly, we propose that Akin scales 
as r 2 : 

Akin oc r 2 . (17) 

Under this scheme, the main cause for evolution is the spin down of the black hole, 
manifested as a decrease of a due to the continuous energy extraction. This way the power 
of the source decreases with cosmic time, giving rise to the observed positive cosmological 
evolution, in the sense that sources were more powerful in the past (the apparent negative 
evolution of XBLs is addressed in 



We additionally assume that the intensive physical variables that describe the 
relativistic jet in BLs, such as the magnetic field, the comoving plasma density, and the 
high energy cutoff of the injected electron energy distribution, all vary over a small intrinsic 
range of values. The combination of these two assumptions suggests a unification scheme 
in which the observed properties of a BL depend mainly on the electron kinetic luminosity 
Ak in of the jet and the angle between the line of sight and the jet axis B. A similar scaling 
relation has been discovered observationally between the size Rbelr of the broad emission 
line region (BELR) and the central (accretion disk) ionizing luminosity L acc of Seyfert 
galaxies and quasars flKaspi et al. 1996) ) using reverberation mapping: 

L acc oc Rbelr- (18) 

This scaling law implies an ionization parameter U (ratio of ionizing photon density to 
electron density) that has a small intrinsic range of values (U ~ 0.1 — 1) for objects that 
extend in luminosity over almost four orders of magnitude ( [Wandel 1997|) . Observational 
arguments, suggestive of a self-similar description of radio-loud AGN, have been presented 
by Rawlings and Saunders (1991), who find that the kinetic luminosity of the kpc-scale 
jet is proportional to the narrow emission line luminosity. This relation holds for objects 
with or without broad emission lines and extends for more than four orders of magnitude 
in luminosity. In a similar study, Celotti, Padovani, & Ghisellini (1997) find that Lbelr 
and the pc-scale kinetic luminosity of the jet for a sample of radio loud AGN are of the 
same order of magnitude. As recent studies of radio loud AGN suggest (e.g. |Koratkar et al 



1998|) , the BELR is probably ionized by thermal radiation from the central engine, rather 



than by the jet radiation. 

We adopt as a working hypothesis a picture in which the ionizing luminosity of the 
accretion disk L acc , the jet electron kinetic luminosity Aki n , and the BELR luminosity 
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Lbelr are closely related, scaling as 

L acc oc L BELR oc A kin oc r^, (19) 

where r* is a dimension that characterizes the jet. We stress here that this scaling 
assumption does not refer to a specific jet model, but is a rather general, model independent 
hypothesis that can be applied to specific jet models. 



5.2. Application to the accelerating jet model 

We apply this self-similar scheme to the accelerating jet model. The question we 
address here is whether we can populate the v s -L s diagram of Figure ||f, under the 
assumptions that the jet electron kinetic luminosity scales with the square of the jet size, 
and that the physical parameters which describe the jet have a small intrinsic range of 
values. We perform calculations using two somewhat different families of values for the 
physical parameters, listed in Table [5]. For each family we run the simulation for a range of 
angles and luminosities, scaling the square of r± and z+ with the electron kinetic luminosity 
Akin- The results of this set of simulations can be seen in Figure || The solid curves 
represent calculated values of L s for a constant electron kinetic luminosity and size as the 
angle varies between 0° and 18°. The dash-dot curves correspond to constant angle as 
the electron kinetic luminosity A kin increases from A to 256A, where A = 10 43 erg s" 1 . This 
increase in luminosity is accompanied by an increase in the jet size: r* and increase by a 
factor of 16 , starting at r* = 2.5 x 10 14 cm and 2* = 5.0 x 10 14 cm. 

The qualitative behavior of both model families is the same. For a given electron 
kinetic luminosity, as the angle increases, the peak frequency v s increases and the peak 
apparent luminosity L s decreases. For a constant angle, as Akin decreases, there is a drastic 
decrease in L s and a milder increase in v s . The two families of models generated by these 
somewhat different physical descriptions cover a significant portion of the observed range 
of values. This behavior is not specific to the above choice for the physical parameters: a 
different set of jet parameters would generate a pattern similar to those obtained by the 
models of Fig. [5]/Table |3[ Modest changes in the physical description merely move the grid 
around the observed parameter space. A small range of values for the physical description of 
the jet, coupled with the 0-A scheme, can populate the v s -L s diagram for BLs. Choosing 
the appropriate set of jet parameters, the 0-A scheme can produce synthetic spectra of 
BL18-BL19 sources. We choose not to do so, since the determination of u s for the two 
sources with \ogv s > 18 in Sambruna et al. (1996) is problematic. 

The conclusion reached previously by Sambruna et al. (1996) and Georganopoulos 
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and Marscher (1996) that the range of v s cannot be explained merely by jet orientation is 
supported by this study. However, a combination of jet orientation and a relatively small 
range in the physical parameters that describe the jet can explain the range of v s under 
the 0-A scheme. A quantitative comparison, including specific assumptions about the 
range of the physical parameters and extensive comparison of synthetic samples to existing 
samples, does not seem worthwhile until the parameters can be more strictly constrained, 
for example through analyses of multiwavelength variability. 

One might naively expect that, for a given angle, since the jet electron kinetic 
luminosity and size scale in a self-similar fashion, the SED would follow this self-similar 
behavior by scaling up or down, following the electron kinetic luminosity and keeping its 
shape and peak frequency the same. An important feature emerging from our simulations 
is that for a constant angle, as the electron kinetic luminosity and the size of the jet 
increase, the peak frequency decreases. This means that, although the jet scales in a 
self-similar manner, this self-similarity is not transfered completely to the SED, since the 
peak frequency changes. The reason for this break-down is that the effects of synchrotron 
losses (cf. eq. [12] and the discussion of §§ |3.1|) on v s do not follow scaling laws that are 



qualitatively similar to those that govern the behavior of L v at v < v s . For example, if 
the values of r* and were doubled while holding B± constant, v s would decrease because 
the electrons would suffer more severe synchrotron losses as they traverse the characteristic 
distance before encountering the region of the jet where B decreases. 



5.3. The negative evolution of XBLs 

The evolution of all the AGN families is positive, in the sense that these objects 
were either more common or more luminous (or both) in the past. A unique and striking 
exception is the negative evolution of XBLs: Morris et al. (1991) showed that the evolution 
of the EMSS XBL sample is negative. It seems that these objects were either less common 
or less luminous in the past. This has been verified recently by Perlman et al. (1996) for 
the ROSAT PSPC sample of XBLs, with (V e /V a ) = 0.331 ± 0.060, where V e is the volume 
of a sphere with radius equal to the distance to each object and V a is the available volume 
within which each object could have been detected in the EMSS. This result is in contrast 
to the value (V e /V a ) = 0.60 ± 0.05 obtained by Stickel et al. (1991) for the 1 Jy RBL 
sample. 

This negative evolution of XBLs can be explained naturally in the framework of the 
G-A scheme, under the assumption commonly invoked for AGN that the electron kinetic 
luminosity Akin of the radio-loud AGN jets has a positive evolution, viz. jets (on average) 
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in AGN had a higher A^n in the past. As we showed in § §5.2] for a given set of physical 
parameters and angle 0, as Akin decreases, the peak synchrotron apparent luminosity L s 
decreases and the peak synchrotron frequency u s increases, essentially shifting the source 
from the region of the RBL-type objects toward that of the XBL-type objects. If we start 
with a sample of sources characterized by a range of A^n and and let these sources 
evolve to lower values of electron kinetic luminosity A^ n , then each source will gradually 
(with cosmological time) shift its characteristics to become more XBL-like. This means 
that the fraction of sources that would be classified as RBLs will gradually decrease with 
time/redshift, resulting in (V e /V a ) > 0.5 in a radio-flux limited sample, while the fraction of 
sources that would be classified as XBLs will gradually decrease and yield (V e /V a ) < 0.5 in 
an X-ray flux-limited sample. 

Our model predicts that the value of (V e /V a ) of an intermediate BL sample, such as 
that of Laurent-Muehleisen (1997), will be confined between the corresponding values of 
(V e /V a ) for the XBL and RBL samples mentioned above, exhibiting almost no evolution 
{(Ve/Va) ~ 0.5), since the population of intermediate BLs on one hand will gain members 
from the RBLs, and will lose members to the XBLs as cosmic time increases. In fact, 
a very recent report ( Bade et al. 19981 ) agrees with this prediction. Bade et al. (1998) 
have studied a new sample of XBLs selected from RASS. They subdivide their sample 
into two halves according to the X-ray to optical flux ratio and find that the extremely 
X-ray dominated subgroup shows negative evolution, while the subgroup with intermediate 
SEDs is compatible with no evolution at all. They also find, however, that the extremely 
X-ray dominated subgroup has higher X-ray luminosity and redshift than the intermediate 
subgroup, an effect that does not seem to agree with our model. 



5.4. Emission line strength under the 0— A scheme 



The discussion in §§ |5.1| implicitly treats BLs as AGNs with a BELR and a central 
isotropic source of thermal ionizing radiation, essentially reducing the differences between 
BLs and FSRQ. The original definition of a BL Lac object includes the criterion of absent 
or weak (equivalent width W\ < 5 A) emission lines QStickel et al. 199Tj ) . This has led 



most workers in the field to assume that either the gaseous environment near the central 
engine of BLs is very poor, or there is no significant source of ionizing radiation, or both. 
Additionally, the conviction of a bimodality separating BLs and FSRQ in terms of their 
emission line properties had been established. The situation, though, is far from clear and 
recent results point toward a continuous transition from BLs to FSRQs in terms of their 
emission line characteristics. Even in the complete sample of Stickel et al. (1991), the 
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authors identified 6 out of 34 RBLs that had emission lines violating the equivalent width 
limit of 5 A. It is also true that the equivalent width is variable and there are objects 
that change classification depending on the time of observation ( [Antonucci 1993j ; |Ulrichj 
Maraschi, & Urry 1997|) . 



The BL PKS 0521-365 flScarpa, Falomo, fc Pian 1995|) , and the prototype object of 



this class, BL Lacertae ( |Vermeulen et al. 1995| ), have shown emission lines with equivalent 



widths W\ > 5 A . More surprisingly, the X-ray selected BL Mkn 501, which has recently 
been detected at TeV energies flQuinn et al. 1996]) , has shown emission lines in the past 
( |l\loIes, Masegosa, fc del Olmo 1987]) . Corbett et al. (1996) argue in the case of BL Lacertae 



that the most plausible way to power the emission lines is through thermal radiation from 
an accretion disk. The required thermal radiation power does not alter significantly the 
total observed flux and the object does not show the "big blue bump" component that 
many radio loud quasars exhibit. 

In a recent study of emission-line luminosities for a sample of HPQs and RBLs, Scarpa 
& Falomo (1997) find that the line luminosity ranges of HPQs and BLs largely overlap, 
although RBLs have smaller equivalent widths. Specifically, in Figure 10 of their paper they 
plot the line luminosity in Mg II versus the continuum luminosity for the sources of their 
sample, as well as the line that corresponds to the equivalent width limit W\ = 5 A. Their 
results clearly demonstrate that for any given continuum luminosity the corresponding 
emission line luminosity does not show any gap separating RBLs and HPQs. If we were to 
increase the arbitrary 5 A limit then we would accept as BLs sources that previously have 
been classified as HPQs. Similarly, lowering the equivalent width limit would shift objects 
from the BL to the HPQ family. It therefore seems plausible that the gaseous environment 
of BLs, although poorer than that of the HPQs, may not be qualitatively different. 

Additional arguments for the presence of ionized gas close to the central engines of BLs 
comes from studies of X-ray and EUV absorption features in XBLs. Madejski et al. (1991) 
found in all objects from a sample of five XBLs an absorption feature at an energy of ~ 650 
eV. Konigl et al. (1995) studied the X-ray bright BL PKS 2155-304 with the Extreme 
Ultraviolet Explorer and found an absorption feature between ~ 75 and ~ 85 A. They 
modeled this feature as absorption resulting from gas clouds at a distance from the central 
engine similar to the typical BELR size, moving with a speed ~ 0.1 c. Similar results have 
been reported recently from EUV observations of the XBL Mkn 421 ( [Kartje et al. 199"T ), 



and X-ray observations of the XBL H1426+428 ( [Sambruna et al. 1997| ). 

In Figure |] we plot the synchrotron peak apparent luminosity L s as a function of the 
synchrotron peak frequency v s for the XBL (diamonds) and RBL (triangles) samples of 
Sambruna et al. (1996). The RBLs with emission lines that violate the 5A equivalent width 



-22 - 



criterion are marked with small filled circles. Note that these sources are preferentially 
low u s , high L s objects. We use the scaling relation flT9|), coupled with the accelerating jet 
model to show that the lack of emission lines in high u s , low L s (BL15-BL17) BLs and 
the existence of emission lines preferentially in low u s , high L s (BL13-BL14) BLs can be 
reconciled with the existence of a scaled-down accretion disk illuminating a BELR. 

In Figure we plot the SED for a range of A kin . We assume that the BELR luminosity 
is proportional to the jet electron kinetic luminosity, and for the purpose of demonstration 
we consider a broad line with a full width at half maximum (FWHM) of 3000 km/s and 
with power Lbelr = 0.3Afcj n . The BELR luminosity is shown by the short horizontal line 
centered at v — 10 15 Hz. As A kin and L BELR decrease, the synchrotron peak frequency u s 
increases. Although the decrease of Kkin is accompanied by a decrease in the emitted power 
vL v at every frequency, the shift of v s toward higher frequencies results in the gradual 
dominance of the synchrotron emission over the line radiation. In Figure ||] we plot as filled 
circles the loci of the model for the range of luminosities used in Figure |7j. The qualitative 
behavior of the model is similar to the trends exhibited by the data: as the model points 
shift to higher values of v s and lower values of L s , the equivalent width W\ decreases, 
crossing over to values smaller than 5A. 

Thus, although the ratio of jet electron kinetic luminosity to BELR luminosity remains 
the same, the equivalent width of the emission lines decreases as the electron kinetic 
luminosity of the source decreases. This is again a manifestation of the breakdown in 



self-similarity described in §§5.2 



6. SUMMARY AND DISCUSSION 

The main points of this work follows. 

1. We verify that an increase in the angle G between the line of sight and the jet 
axis shifts the observed characteristics of the source from the RBL to the XBL class. As 
has been noted, however ( |Sambruna, Maraschi, fc Urry 1996] ; peorganopoulos fc Marscher 



|1996[ ), this shift is not enough to explain the range of observed physical parameters in 
complete BL samples. 

2. We propose a new continuous classification scheme for BLs, based on the synchrotron 
peak frequency v s of the SED of a source. This scheme avoids the observationally induced 
bimodality of the currently used XBL-RBL, or HBL-LBL schemes, and naturally 
accommodates the newly discovered ( [Laurcnt-Muchlciscn 199? ; Laurcnt-Muchlciscn ct alT 



1998| ; [Permian et al. 1998| ) intermediate BLs. As can be seen in the results of Sambruna et 
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al. (1996), given the class of a BL, we can infer upper limits to the synchrotron luminosity 
of the object, the steepness of the spectrum at different frequencies, and the maximum 
concavity index a ox — a x . 

3. Observations suggest that BLs may be characterized by a scaled-down accretion 
disk and BELR, similar to those found in FSRQs. Based on recent observational results, we 
propose a self-similar blazar description, according to which the accretion disk luminosity, 
the BLR luminosity, and the jet kinetic luminosity scale with the square of the scale size 
of the jet [expression ([19])]. This description leaves invariant all the primary physical 
parameters such as the pressure and particle and energy density at the base of the jet. 

We propose the 0-A scheme unification scheme: the electron kinetic luminosity, 
together with the orientation of the jet, determines the observed characteristics of a source. 
A small range of the physical parameters that describe the jet, together with the O-A 
scheme, can reproduce the range of observed BL synchrotron peak frequencies u s , which 
cannot be reproduced solely by orientation effects ( Sambruna, Maraschi, fc Urry 1996 ; 
Georganopoulos fc Marscher 1996| ). 



4. The negative evolution of the XBLs can be explained naturally in the context of this 
self-similarity based scenario. If we start with a random sample of sources, as cosmological 
time passes the electron kinetic luminosity decreases, the peak frequency v s for each source 
increases, and the number of sources that would be classified as XBLs increases. We predict 
that samples of intermediate sources will have a (V e /V a ) value intermediate between those 
of XBL and RBL samples, as new results have indicated. 

5. We show that the existence of emission lines preferably in low v s (BL13-BL14), 
high L s sources, and the lack of lines in higher u s , lower L s objects, is reproduced under 
scaling relation flUf) . As both A^n and Lbelr decrease, v s shifts toward the frequency of 
the emission lines and the non-thermal continuum gradually dominates over the emission 
line luminosity. 

Our discussion of the O-A scheme has been confined to the BLs, although there are 
many indications that a common underlying physical mechanism could describe all blazars. 
Since the luminosity output of FSRQs is dominated in many cases by the 7-ray flux and 
our simulation does not calculate the inverse Compton energy losses, we cannot at present 
apply our model to FSRQs. The extended radio luminosity distribution of FSRQs is similar 
to that of FR II radio galaxies ( |Maraschi fc Rovetti 1994]) , and the jet magnetic field tends 
to be aligned with the jet axis ( |Cawthorne et al. 1993| ). BLs, on the other hand, have an 
extended radio luminosity distribution similar to that of FR I radio galaxies ( |Perlman & 



Stocke 1993| ), and the jet magnetic field tends to be perpendicular to the jet axis flGabuzda 
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et al. 1992). 

Although it may seem as if these differences are enough to separate BLs and FSRQs 
into two distinct families, there are many observational arguments that support the 
possibility of a continuous distribution of sources that belong to a single, continuous family. 
We have already reviewed observations suggesting that a separation of blazars on the basis 
of the emission line properties may not be justified. The extended radio power distribution 
for BLs and FSRQs overlap significantly, and the same holds for their parent populations, 
the FR I and FR II radio galaxies flMaraschi fe Rovetti 19941 ; [Padovani 1992|) . A population 



of FR II radio galaxies with low-excitation narrow-line spectra, low radio power, and more 
relaxed radio morphologies ( |Laing et al. 1994] ) is indicative of a connection between FR 



I and FR II sources. Recent VLBI studies flUawthorne fc Gabuzda 1996| ; Kollgaard et 



al. 1996; |Gabuzda et al. 1994|; [Leppanen, Zensus, fc Diamond 1995|; [Kemball, Diamond 



Sz Pauliny-Toth 1996 ) indicate that the magnetic field orientation does not always follow 
the previously observed bimodality, and FSRQs with perpendicular and BLs with parallel 
magnetic field are not rare. Additionally, recent results ([Lister fc Marscher 1998|) from 
higher frequency (43 GHz) VLBI observations indicate that the magnetic field orientations 
of both BLs and FSRQs in the radio core region are similar. 

There are several BLs in the samples that Sambruna et al. (1996) consider which 



have extended radio characteristics that correspond to FR II radio galaxies (|Kollgaard 



et al. 199^ ; [Murphy, Browne, fc Perley 1993[ [Wrobel fc Lind 1990|) , or have longitudinal 



magnetic fields QGabuzda et al. 1994j) . As can be seen in Figure || all these BLs are low u s 



(BL13-BL14), high L s objects that occupy a region of the parameter space in the v a -L a 
diagram similar to that occupied by the FSRQ sample. It therefore seems justifiable to 
consider that the properties of BLs and FSRQs do not correspond to two distinct physical 
structures, but, rather, form a continuum. 

Since the jets in these sources are relativistic, the angle that the jet forms with the line 
of sight must be one of the parameters that affects their observed properties. The range of 
observed luminosities for objects with similar synchrotron peak frequencies v s has led us to 
propose that a second crucial parameter is the electron kinetic luminosity of the jet. Our 
simulations show that this formulation indeed holds promise for a global understanding 
of blazars. For further progress, detailed, time-dependent simulations that include in an 
accurate fashion the relevant physical processes, such as inverse Compton scattering of 
both jet synchrotron photons and external photons, need to be developed and compared to 
observations of individual sources and to the statistical properties of complete samples. 
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A. GAS DYNAMICS 

We describe here the analytical one-dimensional steady flow described by Blandford & 
Rees (1974). The plasma produced in the vicinity of the central engine is ultrarelativistic 
in the fluid proper frame (mean total to rest mass ratio per particle 7 ^> 1), and remains 
in the ultrarelativistic regime as it flows downstream. The flow is considered adiabatic and 
quasi-spherically symmetric, i.e., that the jet is confined to a cone. The pressure profile is 
parameterized by a power law, 

p(z)=p ^y a (Ai) 

where the z axis is the symmetry axis of the system, z is the distance of the flow 
stagnation point (bulk motion Lorentz factor T = 1) from a fiducial point where the 
pressure mathematically becomes infinite, po is the pressure at zo, and a is the exponent 
that describes how fast the pressure drops (a > 0) as we move away from the central engine 
along the z axis. Note that z = at the fiducial point and not at the central engine. The 
jet is characterized by a constant power, 

A t = T 2 {3wcA, (A2) 

where (5 is the jet velocity in units of c, w is the enthalpy density in the fluid frame and A 
is the cross-section of the jet at a given z. For a relativistic gas w = (4/3)e, where e is the 
comoving energy density of the plasma. 

One way to ensure that the flow is adiabatic is to assume that a significant portion 
of the energy is carried by relativistic protons that have essentially no radiative losses. 
Alternatively, we confine our study to EED with s > 2, which corresponds to most of the 
energy being near the low energy cutoff of the EED. Since the low energy electrons have 
small radiative losses, such a flow is approximately adiabatic. We use the first description: 

e = e e i + e p , e e i = ke, e p — (1 — k)e, < k < 1, (A3) 

where 

eei = / lN^~ s dl = - (A4) 

is the electron energy density and e p is the proton energy density in the comoving frame. 
Assuming equipartition between the electron and proton energy density, we set k = 0.5 
(the actual requirement is k < 0.5). This ensures that half of the energy is in protons that 
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do not radiate, and therefore our adiabatic approximation is reasonable. Therefore, the 
electron kinetic luminosity Akin is 



A km = kA t = A t /2, 

The cross-section of the jet is minimized at the sonic point (I\ 
characterized by a pressure 



(A5) 



a cross-sectional radius r\ 



and a distance z± from the fiducial point 

'9\ 1/a 



Z± — Zq 



3/2), which is 

(A6) 

(A7) 
(A8) 



The cross-sectional radius of the jet as a function of z is given by 



3e 

Z \ 2 

/'( :) = ; . I — 

- z* ■ 



z 



2c 



3(-] -2 

z± 



and the bulk Lorentz factor of the flow V is given by 



r(*) = r, 



(A9) 



(A10) 



The exponent e determines how fast the jet opens and accelerates and is determined solely 
by the steepness of the external pressure gradient, e = a/4. For z ^> z+, r oc z e . 



B. THOMSON LOSSES 

Our treatment of electron losses does not consider the Thomson losses due to the 
radiation from an accretion disk. If z x is the distance of the fiducial point used in the flow 
description, the distance of the base of the jet from the central engine is Zbase — z\ + z±. 
The synchrotron losses will be greater than the Thomson losses if the accretion disk photon 
energy density in the comoving frame at the base of the jet is greater than the comoving 
magnetic field energy density. 
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As given in Dermer & Schlickeiser (1994), the comoving photon energy density u acc at 
distance z from the central engine due to a point-like accretion disk is 

L acc 

where L acc is the accretion disk luminosity, and I\, (3* describe the jet flow at the base of 
the jet. Requiring u acc < B+/(8ir) we obtain 

z b ase > z min = (2/3) 1 / 2 10 17 4{ 2 J B- 1 r- 1 (l + cm, (B2) 

where L 44 is the accretion disk luminosity in units of 10 44 erg s _1 . In this work we assume 



that the fiducial point is adequately far from the central engine so that relation |B2 
satisfied. 



is 
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Fig. 1. — The model spectral energy distribution (SED), for values of the parameters as 
given in Table [j], except for the single parameter allowed to vary in each panel. Here, L u is 
the apparent luminosity per unit frequency 
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Fig. 2.— The (bottom panel) SED for values of the parameters given in Table [l], for 
various observing angles G between the jet axis and the line of sight. Peak luminosity L s 
(middle panel) and optical to X-ray spectral index a ox (top panel) versus observed peak 
synchrotron frequency v s . The data points correspond to samples of RBLs (triangles) and 
XBLs (diamonds) flSambruna, Maraschi, fc Urry 1996| ), and the curves correspond to the 
model. 
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Fig. 3. — The spectrum of PKS 2155-304 during the multiwavelength campaign of November 
1991. The data were obtained on November 14 during an 18-hour period. ( [Edelson et al. 



1995|) . The solid curve corresponds to the model fit produced using the parameters listed in 
Table Q 
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Fig. 4. — The (a) radio to X-ray (a rx ), (b) radio to optical (a ro ), (c) optical to X-ray (a ox ) 
broadband spectral indices; (d) the X-ray concavity index a ox — a x ; (e) the bolometric 
synchrotron apparent luminosity Lb', and (f) the synchrotron peak luminosity L s , as a 
function of the observed synchrotron peak frequency v s of the synchrotron power spectrum 
for the XBL (diamonds) and RBL (triangles) samples of Sambruna et al. (1996). The 
arrow indicate an RBL with only an upper limit to the peak frequency v s . The dashed lines 
separate the populated from the unpopulated areas of the diagrams. 
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Fig. 5. — The synchrotron apparent peak luminosity L s as a function of the observed 
synchrotron peak frequency u s of the SED for the XBL (diamonds) and RBL (triangles) 
samples of Sambruna et al. (1996). The two superposed grids are the results of the 
accelerating jet model coupled with the scaling assumption of equation ( |T9"D for two different 
physical descriptions of the jet, as given in Table [3[ The solid curves represent constant 
electron kinetic luminosity Akin, while the dash-dot curves correspond to constant angle G 
between the jet axis and the line of sight. 
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Fig. 6. — The synchrotron apparent peak luminosity L s as a function of the observed 
synchrotron peak frequency v s for the XBL (diamonds) and RBL (triangles) samples of 
Sambruna et al. (1996). Small filled circles mark RBLs with emission lines that violate 
the 5A equivalent-width criterion. The bigger filled circles, connected with a dotted line, 
represent the model points of Figure [?|. The most luminous model source is assigned an 
arbitrary equivalent width W\ = 7 A , and the variation of the equivalent width W\ is shown 
as Akin and Lbelr decrease. 
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Fig. 7. — The SED for a range of jet electron kinetic luminosities. The BELR luminosity, 
shown by the short horizontal line is proportional to the jet electron kinetic luminosity (see 
§§ p.4| ). The parameters that describe the jet are the same as the parameters describing the 
family A jets in Table |3|, with = 3°. 
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Fig. 8. — Peak synchrotron apparent luminosity vs. observed peak frequency for the two 
BL samples (see Fig. ||) and the small complete sample of FSRQ (filled squares) from the 
S5 survey (Brunner et al. 1994), used by Sambruna et al. (1996). Small filled circle mark: 
(a) RBLs with magnetic field parallel to the jet (0735+178, 1308+326, and 2007+777 from 
Gabuzda et al. 1994); (b) RBLs with FR II like extended radio characteristics (0235+164, 
1308+326, 1538+149, and 1823+568 from Murphy et al. 1993, 0954+658, 1749+701, and 
2007+777 from Kollgaard et al. 1992, and 1807+698 from Wrobel et al. (1990); (c) RBLs 
with emission lines that violate the 5A equivalent width criterion (0235+164, 0537-441, 
0851+202, 1308+326, and 1749+096 from Stickel et al. (1991), and BL Lacertae from 
Vermeulen et al. 1995). 



-41 - 



Table 1. Model parameters. 



Parameter Value 



r* 10 15 cm 

z± 2.0 x 10 15 cm 

r* 2.0 

e 0.3 

0.5 G 

m 1.0 

A fel „ 10 45 crgs/s 

7ma:c 2.0 x io 5 

s 2.0 

e io° 



Table 2. Model parameters for PKS 2155-304. 



Parameter 


Value 


r± 


1.1 x 10 15 cm 




1.5 x 10 15 cm 


r* 


1.5 


e 


0.3 




0.11 G 




8.63 x 10 45 ergs/s 


'Jmax 


6.5 x 10 5 


S 


1.7 


e 


10° 
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Table 3. Values of parameters used in simulations under the 0-A scheme (cf. Fig. |5|). 



Parameter 


Family A 


Family B 


r* 


rf 2 2 r, 2 4 r 


rf 2 2 r, 2 4 r 


z* 


z*2 2 z, 2 4 z 


z^2 2 z, 2 4 z 


r* 


1.5 


2.0 


e 


0.3 


0.4 


B* 


0.4 G 


0.4 G 




A^ 2 4 A, 2 8 A 


Af 2 4 A, 2 8 A 


^fmax 


4 x 10 5 


8 x 10 5 


S 


2.0 


1.7 



a r = 2.5 x 10 14 cm 
h z = 5.0 x 10 14 cm 
C A = 10 43 ergs/s 



